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Abstract. A review of the basic principles of population synthesis for binary stars is presented. 
We discuss the break-up of low and intermediate mass close binaries over different evolutionary 
scenarios and, as an example, briefly consider results of the population synthesis for SN la. 



1. INTRODUCTION 

Population synthesis for binary stars is a convolution of the statistical data on initial 
parameters and birthrates of binaries with the evolutionary scenarios for them. The main 
goals of population synthesis are: 

• to understand the descent of binaries of different types from unevolved main- 
sequence binaries which differ only in the masses of components A/i,A/2, and their 
separations a and to understand the links between the objects of different types; 

• to estimate the incidence of binaries of different types and occurrence rates of certain 
events in them (for instance, explosions of Novae or Supemovae); 

• to find the distributions of binaries of different types over parameters like masses of 
components, orbital periods, luminosity of components, spatial velocities; 

• to model 'observed" samples of stars applying selection effects. 

Below, we review the basic principles of population synthesis, discuss the break-up of 
low and intermediate mass close binaries over different evolutionary scenarios and, as 
an example, consider results of population synthesis studies for Supemovae la (SNIa). 
We define low and intermediate mass binaries as objects in which components end their 
evolution as white dwarfs. 



2. EVOLUTIONARY SCENARIOS FOR BINARY STARS 

The main concept of the population synthesis is that of an "evolutionary scenario", the 
sequence of transformations of a binary system with given initial (Mio,M2o,«o) that it 
can experience in Hubble time (the term is coined by van den Heuvel and Heise (?]). 

For construction of evolutionary scenarios all binaries are separated into "close" and 
"wide" ones. In close binaries components interact between themselves, the principal 
form of interaction is mass exchange through Roche lobe overflow. The possibility 
for interaction of components in a given binary is defined by critical radius Re, a 
boundary within which the atmosphere of the evolving star can expand while star 
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FIGURE 1. Descendants of components of close binaries depending on the radius of star at the instant 
of RLOF. Close in evolutionary sense stars fill Roche lobes before expanding to /?max- For wide binaries 
the products of evolution are the same as for single stars. 



remains a distinctly separate object. If R^, > Rcr, stellar matter starts to flow out. Figured] 
sketches ZAMS- and maximum radii of stars and shows descendants of components of 
binaries depending on their mass and radius at the instant of RLOF (for solar chemical 
composition stars). The upper mass limit for progenitors of He white dwarfs (WD) is 
certain to several tenth of Mq. The boundary between progenitors of WD and neutron 
stars (NS) is known to about IMq (it may be different for single stars and components of 
close binaries). The least certain is the lower cut-off of the masses of progenitors of black 
holes (BH). It is expected that between 20 and 50Mq, components of close binaries 
may produce NS or BH, depending on such parameters as stellar winds, rotation, and 
magnetic fields the limit of 40 M© is sketched in Fig. 1 provisionally. 

The data necessary for construction of evolutionary scenarios is provided by stellar 
evolution computations. The main data are: the rates of mass loss through RLOF or 
stellar wind (SW) for stars in different evolutionary stages and the nature of the products 
of evolution; time scales of evolutionary stages; initial-final mass relations; response of 
stars to accretion; transformations of separations due to mass exchange/loss and angular 
momentum loss, especially in the common envelope stages. 

Evolutionary scenarios for binaries with certain initial parameters may be convolved 
with star formation rate, binarity rate, distributions of binaries over Mio, ao, qo = 
M20/M10 to get birthrates of stars of different types. Convolved with lifetimes this gives 
number distribution of stars of different kinds at any epoch of the Galactic history, their 
distributions over basic parameters, and occurrence rates of different events, e. g., SN. 
Applying selection effects one gets "observed samples" of different binaries. The first 
numerical population synthesis studies were carried out independently by Kornilov and 
Lipunov [5] for NS, Dewey and Cordes [6] for radiopulsars, Politano |iZJ for cataclysmic 
variables, Lipunov and Postnov ||^ for compact magnetized stars. 



Binarity rate for stars is not known for sure, it was studied for separate groups of stars 
only. Studies of selection effects for spectroscopic and visual binaries \^ suggest 
that binarity rate may be close to 100%. 

The studies of spectroscopic and visual binaries suggest that (i) IMF for 

primary components of binaries may be well approximated by power law dn / dM 
M^^-^; (ii) distribution of binaries is flat in log a. This allows to estimate that among 
binaries ~ 40% are close and ~ 96% are low and intermediate mass ones. 



2.1. Common envelopes 

The crucial problem that plagues all population synthesis studies is uncertainty about 
mass and angular momentum loss (AML) from the system and especially about common 
envelopes (CE). It is assumed that CE arise when accretor in a close binary is unable 
to accrete all matter supplied by the donor [11]. The usual assumption is that the 
envelope of the donor engulfs both components; the latter spiral in, losing energy. 
This results in reduction of the separation of components and may end either in the 
merger of components or in the dissipation of CE. This is understood qualitatively, but 
neither analytical theory nor convincing enough computations of the processes inside 
common envelopes exist. Therefore, usual approach is based on comparison of the 
energy released in rapprochement of components and binding energy of CE by applying 
(most commonly) an equation 

M-AMi-Mr) f Mrm Mim\ 

T?ce ^ ^ h (1) 



where indexes i and / refer to initial and final masses and separations, m is the mass 
of accretor, rjce is the efficiency of deposition of energy into common envelope, A is a 
function of stellar radius that, in essence, describes how well donor envelope is bound 
to the core (it varies ~ 3 — 4 times along RG and AGB and depends on the definition of 
"core" boundary and possible release of thermal energy u4, 15, T^]. If whole ensemble 



of close binaries is considered and Eq. ([T]) is applied to any system, irrespective to the 
evolutionary state of components, satisfactory agreement with observations is usually 
achieved if T^ce-^ :^ I —2. 

However, Nelemans and co-authors iflTj] , for instance, have shown that it is impossible 
to reconstruct evolution of observed close binary helium WD back to main-sequence 
binaries if Eq. ([T]) is applied for the stage of unstable mass transfer between components 
of comparable mass. Instead, it was suggested in UtIi to use for this stage an equation 
based on the angular momentum (7) balance: 

-r = r , K (2) 



J \Md+Ma 

where AM is the amount of matter lost by the system. Nelemans and Tout lUsIl have 
shown that applying Eq.( |2l) with 7 ~ 1.5 it is possible to explain the origin of all 
known binary WD and pre-CV. Equations ([T]) and @ give similar results if masses of 
components strongly differ [18J. 



In the case of super-Eddington accretion in the systems with NS or BH accretors 
formation of common envelopes, perhaps, may be avoided. The energy of the accreted 
matter released close to the compact object (at rc) may be sufficient to expel the "exces- 
sive" matter from the Roche-lobe surface around the compact object (at tl.c) if accre- 
tion rate M < MEdd''L,c/''c (see, e. g., [19]). It is also suggested that common envelopes 
may be avoided by driving excess of matter by optically-thick winds from the surface of 
hydrogen-burning accreting white dwarfs [20]. Note, that the model of radiatively driven 
winds still lacks elaborated theoretical justification, but it is supported by the presence 
of outflows in many systems with compact components. 

If CE formation is avoided as described above, it is usually assumed that the matter 
leaves the system with specific angular momentum of the accreting component. This 
assumption proved to be useful, for instance, for the models of formation of millisecond 
radiopulsars 1E1II22I1 or ultracompact low-mass X-ray binaries with WD donors [ Esj] . 



2.2. Evolutionary paths of low and intermediate mass close binaries 

The character of mass-exchange in a binary depends on the mass ratio of components 
and the structure of the envelope of Roche-lobe filling star [24, 25]. For stars with 
radiative envelopes, to the first approximation, mass exchange is stable if ^ ^ 1.2; for 
1.2 ^ ^ ^ 2 it proceeds in the thermal time scale of the donor; for ^ ^ 2 it proceeds 
in the dynamical time scale. Mass loss occurs in dynamical time scale if the donor has 
deep convective envelope. In practice this means that overwhelming majority (~ 90 %) 
of close binaries pass through one to four stages of common envelope. 

A run of population synthesis code (e. g., ll26t] ') for ~ 1.5 x 10^ close binaries with 
0.8 < M\q/Mq < 11.4 produces ~ 600 different scenarios of evolution. In Fig.|2lwe 
present a kind of summary of scenarios for systems that have first unstable RLOF. For 
description, we denote evolutionary routes in Fig.|2las "path N". 

In the tightest systems primaries overflow Roche lobe when they are still in MS. 
The least massive of these systems (Mio ^ I.5M0) evolve under the influence of 
AML via magnetically coupled stellar wind (MSW) and gravitational waves radiation 
(GWR). An outcome of evolution may be formation of a WUMa system that merges 
into rapidly rotating single star (path 1). About 2% of these systems do not finish their 
evolution in Hubble time and remain MS-stars or giants. Majority of merger products 
end their evolution as WD (path 2-3). Peculiarly, some merger products accumulate 
mass ^ IOM0 and explode as SN II (path 4). 

In wider systems donors overflow Roche lobes when they have He- or CO-cores. 
About one third of binaries avoid merger in CE and form pairs consisting of a He- or 
CO/ONe-WD and a low-mass companion that does not evolve past MS or subgiant stage. 
The systems remain wide enough to avoid merger due to AML via GWR and/or MSW 
(path 6). Systems that merge in CE typically finish their evolution as WD (path 7-8). 

If a binary with a He-core donor does not merge in CE, it first forms a He-star-i-MS-star 
pair (path 9), that may be observed as an sdB star with MS companion. When He-star 
finishes its evolution, a pair harboring a CO white dwarf and a MS-star appears (path 10). 
Similar systems are produced by binaries where RLOF occurs at AGE (path 11). 
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FIGURE 2. Flowchart of the evolution for close binaries that experience unstable first RLOF. The 
numbers indicate "path" discussed in the text. For some important channels the fraction of all close 
binaries passing through them is indicated. 



The tightest of the latter systems (a ~ several Rq) that have low-mass ( ^ I.SMq) MS- 
companions to WD evolve under AML and may form cataclysmic variables (path 12). If 
accreted hydrogen burns at the surface of WD stably, the system may appear as supersoft 
X-ray source (SSS in the plot). The outcome of the evolution of CV is not completely 
clear. The donor may be disrupted when its mass decreases below several hundredth of 
Mq ll27h : accretor may accumulate enough mass to explode as SNIa (see below); WD 
may experience an accretion induced collapse into a NS and a low-mass X-ray binary 
may be formed. Actually, for a considerable fraction of CV evolution is, most probably, 
frozen when M becomes low. 

The second common envelope may form when MS-companion to WD overfills 
its Roche-lobe (path 13). If the donor has a He-core and system does not merge, 
an sdB-i-WD system may emerge (path 14). If separation of components is suffi- 
ciently small, AML via GWR may bring He-star into contact while He is still bum- 
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FIGURE 3. Flowchart of the evolution for close binaries that experience stable first RLOF. 



ing in its core. If A/hp/Mwh \.1 stable mass exchange is possible with typical 
M ~ 10"^ M© yr^^pS]. It was suggested that accretion of He onto CO WD at such 
M leads to SN, when accreted He-layer detonates and initiates detonation of C (edge-lit 
detonation [29]). However, it seems more likely that the lifting effect of rotation reduces 
the power of He-ignition and a He-nova is produced instead [30] . If the system survives 
SN or nova explosion, an AM CVn-type system may be formed. 

Some 12% of close binaries merge in the second CE and form single WD (path 15- 
16). This may create certain excess of massive WD. If the system avoids merger, a pair 
of WD is formed (path 17). The closest of them may be brought into contact by AML 
via GWR. If in a CO-i-He white dwarf pair conditions for stable mass exchange are 
fulfilled, an AM CVn system forms (path 18), see for details [31, 32, 33]. About 5% of 
close binary WD merge (path 19); it is not yet clear how the merger proceeds and it is 
possible that for He-i-He or CO-i-He pairs a helium star is an intermediate stage ll34l] . 

A fraction of close binaries produces merging CO WD pairs ("double-degenerates, 
DD") with total mass greater than Chandrasekhar mass Mch(path 20) that may explode 
as SNIa id QflQ^J- They are probably the most promising candidates for SNIa (see 
below). 

Figure |3] shows typical evolutionary scenarios for the systems that stably exchange 
mass in the first RLOF. It applies to ~ 10% of all close binaries. In the first RLOF 
(path 1) binaries may be identified with Algol-type systems. Algol phase ends by 
formation of WD-i-MS-star system directly (path 3) or with an interfering phase of He- 



star+MS-star (similarly to (path 9-10) in Fig.El). In rare cases WD+MS remain "frozen" 
in this state (path 3). RLOF by the former secondary usually results in CE, since ^ ^ 1 
(path 4). If the system did not merge in CE, the donor had a nondegenerate He-core 
and system is close enough, a semidetached WD+He-star system may arise due to the 
AML via GWR and in such a system a SN or He-nova may explode (path 5). In this case 
formation of an AM CVn star also is possible, like in (path 14) in Fig. El For the systems 
that avoid merger in CE the following outcomes of evolution become possible: formation 
of a non-merging in a Hubble time pair of WD (path 6); formation of a DD that merges 
due to AML via GWR and possibly produces SNIa or a massive WD (path 7). Merger 
in CE also produces a massive WD (path 8). 

A special case of evolution, that involves binaries with both low/intermediate and high 
mass components, leads to formation of recycled pulsars. First, in the system with initial 
qo ~ 10 the primary becomes a NS. If the system is not disrupted by SN explosion, 
low-mass secondary after a stage of stable RLOF (when the system is a low-mass X- 
ray source) becomes a He WD. In another scenario qo is larger and, before forming a 
WD companion to NS, the system passes through CE. In both cases NS are recycled by 
accretion. 

Another peculiar example of evolution involves binaries with ~ 1 and mass of 
components ^ (6 — 1)Mq. If the first mass exchange is stable, initially less massive 
component accumulates enough mass to explode as SN II. In this case massive white 
dwarf (former primary) gets a young neutron star companion (see ll37L l38ll for details). 



3. CANDIDATE SUPERNOVAE lA 

Figure IH summarizes the main evolutionary scenarios that are currently considered 
to result in formation of candidate SNIa systems. For completeness, the channel of 
symbiotic stars is shown, though, they are not "close binaries" strictly. A detailed review 
of the evolution to SN la is recently published in ll39ll . Here we give a brief outline. 

DD-merger scenario (A) operates both in young (~ 10^) yr and old (~ 10^^) yr pop- 
ulations [26, 39] and may provide for our Galaxy the rate of SNIa (Vja) comparable to 
the one inferred from observations |40]: (4±2) x 10^-^ yr^^ The DD- scenario recently 
got strong support both from the theory and observations. Population synthesis ll4lll pre- 
dicted that it may be necessary to study up to ~ 1000 field WD with V ^ 16 17 for 
finding a proper merger candidate. This was done within SPY-project and resulted in 
discovery of the first super-Chandrasekhar DD {Mt = (1.46±O.12)M0] that will merge 
in 4 Gyr and two DD with mass close to Mch that will merge in one to two Hubble times 
p5]. This means that formation of DD candidate SNIa is quite possible. On the "theoret- 
ical" side, SPH calculations of DD-mergers with account for nuclear burning [34] have 
shown that in the core-disc structure that forms, in the region of interaction the tempera- 
ture rises rapidly, lifts degeneracy, expands the matter and thus quenches thermonuclear 
flash. Flame does not propagate to the center, contrary to previously expected (e. g., 
[|43J). Stable accretion from the disk is possible. On the other hand, deposition of the an- 
gular momentum into the core spins-up rotation of WD, causes deformation of WD and 
AML by distorted configuration via GWR and finally makes possible close-to-center 
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FIGURE 4. Main evolutionary channels resulting in candidate SN la systems. 



explosive ignition of carbon 1144114! 

Scenario B operates both in young and old populations, but gives a minor contribution 
to Via since typical total masses of the systems are well below Mch- 

"Single-degenerate" (SD) scenario (C) is often considered as the most promising 
one. Population synthesis estimates lend to it from ~ 10% [46] to ~ 50% [47] of the 
observed Via. The uncertainty in the estimates of Via from this channel stems from still 
poorly constrained efficiency of mass accumulation in the process of accretion that 
is accompanied by hydrogen and helium burning flashes I 48l l49l Isoll . The estimates 
become more favorable if mass transfer may be stabilized by mass and momentum loss 
from the system. "Observational" problem with this scenario is the absence of hydrogen 
in the spectra of SNIa, while it is expected that ~ 0.15 Mq of H-rich matter may be 
stripped from the companion by the SN shell [51]. Recently discovered SNIa 200 lie 
and similar 1997cy with signatures of H in the spectra [52] may belong to the so-called 
SN 1.5 type or occur in symbiotic systems IsbIi . Also, under favorable for this scenario 
conditions, it seems to overproduce supersoft X-ray sources ll4^ . However, SD-scenario 



got recently support from discovery of a possible companion to Tycho SN 

Scenario D may operate in populations where star formation have ceased no more 
than ~ 1 Gyr ago and produce SN- scale events at the rates that are comparable with 
the inferred Galactic SNla rate. But "by construction" of the model, the most rapidly 
moving products of explosions have to be He and Ni; this is not observed. The spectra 
produced by ELD are not compatible with observations of the overwhelming majority of 
SNla 1I55I] . On the theoretical side, it is possible that lifting effect of rotation that reduces 
effective gravity and degeneracy in the helium layer may prevent detonation [30]. 

Scenario E is the only way to produce SN la in a wide system, via accumulation of 
a He layer for ELD or of Mch by accretion of stellar wind matter in a symbiotic binary 
i56]. However, its efficiency is probably low (VTa~ 10^^ yr^), because of typically low 
mass of WD and low M in symbiotic systems iIstIi . 



4. CONCLUSION 

For the conclusion, let us note that, despite evolution of close binary stars is, at least 
qualitatively, understood and for certain systems population synthesis provides even 
quantitative agreement with observations, there are some important problems to solve: 

• mass and angular momentum loss from the systems (including stellar winds), evolution 
in common envelopes; 

• processes during merger of stars, from main-sequence stars through relativistic objects; 
evolution of merger products; 

• the role of rotation, angular momentum deposition during accretion and associated 
instabilities; 

• mass limits between white dwarfs/neutron stars/black holes. 

The author acknowledges financial support from the LOG of the conference and from 
Russian Foundation for Basic Research (grant no. 03-02-16254). 
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